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ABSTRACT 

In  large  electronic  systems  where  the  thermal  loads  are  distributed  in  a  pattern,  coolant  flow  requirements  also  need  to 
follow  similar  pattern.  Amilable  literature  mainly  focuses  on  uniform  flow  distribution  through  manifolds.  However,  for 
rack  mounted  electronics  where  heat  loads  are  in  a  repetitive  pattern,  flow  distribution  also  follows  similar  pattern.  This 
paper  presents  a  methodical  approach  in  providing  the  pattern  flow  distribution  through  manifold  system.  Based  on 
principle  of  pressure  driving  the  flow,  the  desired  flow  is  achieved  by  introduction  ofrequired  pressure  drop  by  means  of 
suitable  orijices.  Three  types  of  patterned  flows  are  presented  in  this  paper.  The  average  percentage  deviations  between 
the  required  flow,  CFD  and  experimental  results  are  compared  and  found  to  be  in  close  agreement  with  each  other.  Based 
on  the  results,  it  is  observed  that  the  proposed  methodology  provides  an  effective  and  simplified  solution  that  can  be 
adopted  in  electronics  cooling  in  case  of  dissimilar  heat  load  characteristics  with  repeated  patterns. 
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NOMENCLATURE 

a:  Cross  sectional  area  of  orifice,  m2 

A:  Cross  sectional  area  of  branch  channel,  m2 

Cd:  Orifice  coefficient  of  discharge 

Cp:  Specific  heat,  kJ/KgK 

d  :  Orifice  diameter,  m 

D:  Branch  channel  diameter,  m 

H:  Heat  load,  W 

m  :  Mass  flow  rate,  kg/s 

Q:  Water  flow  rate,  m3/s 

P:  Volume  flow  rate,  m3/s 

p:  Density,  kg/m3 
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Ap  :  Pressure  drop  across  orifice,  Pa 

At:  Temperature  gradient,  °C 

1.  INTRODUCTION 

Electronic  devices  are  penetrating  all  walks  of  life  and  their  performance  requirements  are  ever  increasing  to  make  them 
more  powerful,  miniaturised  and  reliable.  Since  the  hot  spot  temperatures  affect  component  reliability,  thermal 
management  of  these  devices  by  way  of  maintaining  them  below  the  limiting  temperature  is  an  important  design  focus. 
Different  cooling  methods  are  adopted  to  reduce  hot  spot  temperatures  by  means  of  conduction,  convection  and  radiation 
or  any  combination  of  these  methods.  Cooling  methods  employed  in  electronics  cooling  can  be  air  or  liquid  cooling,  direct 
or  indirect  cooling,  single  or  phase-change  type  [1].  When  the  thermal  loads  are  high,  liquid  cooling  is  resorted  to.  One  of 
the  variants  of  these  cooling  methods  is  indirect  single-phase  liquid  cooling.  In  this,  the  liquid  coolant  circulates  through 
the  heat  sink  and  carries  away  the  heat.  When  numbers  of  electronics  are  more,  coolant  distribution  is  done  through 
manifolds. 

A  manifold  is  a  flow  channel  in  which  fluid  enters  and  leaves  through  the  pores  of  the  channel.  When  the  fluid 
branches  out  through  the  pores  laterally,  it  is  termed  as  distribution  manifold.  In  this,  the  fluid  stream  is  decelerated  in  the 
flow  direction.  When  the  fluid  is  added  from  branch  points  to  the  flow  channel,  in  the  flow  direction,  it  is  called  collection 
manifold.  In  this,  fluid  accelerates  in  the  flow  direction.  Main  flow  channel  and  branch  points  are  also  known  as  header 
and  channels  respectively.  When  both  distribution  and  collection  manifolds  are  connected  and  fluid  direction  is  same  in 
both,  then  it  is  called  parallel  or  Z-configuration  flow  (Figure  la).  When  the  flow  direction  is  opposite  in  the  distribution 
and  collection  manifolds,  it  is  called  reverse  or  U-configuration  flow  (Figurelb). 


Figure  la:  Z-Type  Manifold  Configuration.  Figure  lb:  U-Type  Manifold  Configuration. 

The  flow  configuration  inside  the  manifolds  is  complex  in  nature.  It  is  a  function  of  Reynolds  number,  ratio  of 
sum  of  all  channel  area  to  header  area,  number  of  channels,  pitch  of  channels  and  ratio  of  length  to  diameter  of  the 
manifolds. 

Formulation  of  flow  distribution  inside  the  manifolds  in  the  literature  is  based  on  two  important  assumptions:  (i) 
flow  uniformity  and  (ii)  uniform  pitch  of  the  branch  channels.  Several  discrete  and  analytical  models  framed  by  researchers 
are  on  these  two  assumptions.  These  assumptions  reduce  the  complexity  of  the  governing  equations. 

Flow  distribution  in  manifolds  was  addressed  by  Acrivos  in  1959  [2].  Dutta  &  Majumdar,  Bajura  &  Jones, 
Bassiuony  &  Martin  [3-6]  developed  analytical  models  which  were  validated  experimentally  but  for  manifolds  of  short 
length  by  neglecting  frictional  effects.  J.  Wang  [7-11]  developed  analytical  models  considering  both  momentum  and 


Impact  Factor  (JCC):  8.8746 


SCOPUS  Indexed  Journal 


NAAS  Rating:  3.11 


Numerical  and  Experimental  Imestigations  ofPatterned  Flow  977 

Distribution  through  Manijolds  in  Electronics  Cooling 

friction  factors  along  with  the  flow  inAuencing  factors  in  U-type  and  Z-type  flow  configurations. 

Several  other  researchers  considered  geometrical  modihcations  to  achieve  the  flow  uniformity  in  the  channels. 
These  methods  include  structural  modification  of  headers.  Studies  using  triangular,  rectangular,  trapezoidal  headers 
showed  that  triangular  headers  are  good  in  flow  distribution  compared  to  rectangular  and  trapezoidal  headers  [12].  Also, 
increasing  or  decreasing  the  area  by  tapering  headers  [13,14],  introduction  of  baffle  plates,  changing  the  configurations  of 
increasers  in  the  headers[15]  also  increase  flow  uniformity. 

However,  geometrical  modifications  always  might  not  be  feasible  in  all  applications  in  terms  of  space  constraints, 
and  manufacturing  criticalities.  A  quasi-analytical  method  developed  by  Eprahim  M.  Sparrow  is  very  simple,  yet  having 
potential  application  in  achieving  flow  uniformity  [16].  The  research  by  Zhang  et.  al.[17]  describes  flow  uniformity  in  the 
dividing  manifold  system  with  parallel  pipe  arrays(DMS-PPA)  for  the  anti-parallel  flow  direction. 

All  the  above  studies  reveal  that  performance  of  the  manifold  distribution  has  been  quantified  by  the  degree  of 
flow  uniformity  it  provides.  However,  there  are  situations  where  electronic  units  mounted  in  a  rack  as  shown  in  Figure.2, 
are  of  similar  dimensions,  but  have  patterned  heat  load  distributions.  If  the  flow  distributions  are  designed  considering  the 
highest  heat  load,  it  increases  overall  power,  material  and  space  requirements.  In  such  a  situation,  instead  of  a  uniform 
flow,  a  patterned  flow  distribution  can  provide  an  optimal  solution.  To  the  author’s  knowledge,  very  limited  literature  is 
available  that  addresses  the  desired  or  patterned  flow  distribution. 


Figure  2:  Electronics  Mounted  Inside  the  Rack. 


Formulation  of  analytical  equations  and  finding  solutions  to  this  kind  of  flow  distribution  is  quite  complex. 
Numerical  simulations  to  study  3D  flow  phenomenon  inside  the  manifolds  can  give  more  realistic  approach.  However, 
numerical  simulations  are  specific  to  each  configuration  and  arriving  at  a  suitable  flow  requirement  by  a  unified  approach 
rather  than  trial  and  error  approach  is  the  requirement.  This  paper  provides  a  methodology  of  arriving  at  pattern  flow 
distribution  using  numerical  simulations,  followed  by  experimental  validation.  This  methodology  is  explained  by  using  the 
problem  cases  as  mentioned  in  the  subsequent  section. 

2.  PROBLEM  DEFINITION 

Since  the  objective  is  to  meet  the  flow  requirements,  to  reduce  the  complexity  of  numerical  modelling,  each  cold  plate  is 
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modelled  as  a  flow  channel.  Flow  requirements  to  cool  the  electronics  are  calculated  using  the  formula 

H=mCpAt  (1) 

Where,  Cp  is  a  constant  at  a  particular  temperature.  For  a  defined  temperature  gradient,  mass  flow  rate  of  the 
coolant  is  proportional  to  the  heat  load. 

maH  (2) 

mass  flow  rate  can  also  be  expressed  as 


m  =  pV 

Hence, 

VaR 

Equation  (4)  shows  that,  the  volume  flow  requirements  are  proportional  to  the  heat  load  and  form  a  set  of  repeated 
pattern  for  patterned  heat  loads.  The  calculated  flow  requirements  are  shown  in  Figure3a,  Figure3b  and  Figure3c 
respectively.  All  these  requirements  are  considered  for  the  same  configuration  of  overall  manifold  length  and  pitch  of 
branch  channels. 


(3) 

(4) 


For  each  channel  shown  in  x-axis,  volume  flow  rate  is  shown  as  litres  per  minute  (lpm)  on  the  y-axis.  For  all  the 
cases,  a  total  of  1  lpm  flow  rate  is  considered.  Water  is  considered  as  coolant.  All  the  three  configurations  fall  in  the  regime 
of  laminar  flow  with  Reynold’s  number  <2100. 


Figure  3a:  Pattern-1. 


Channel  No. 

Figure  3b:  Pattern-2. 


Channel  No. 
Figure  3c:  Pattern-3. 
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The  fluid  flow  pattern  across  each  branch  channel  of  the  manifolds  is  governed  by  the  pressure  gradient/drop  across  that 
channel,  inside  the  manifold  system.  Flow  across  channels  is  proportional  to  the  pressure  drop  across  the  channels.  Hence 
by  introducing  appropriate  resistances  in  the  form  of  orifices,  inside  the  flow  channels,  the  desired  pattern  flow  can  be 
attained.  The  sequence  of  proposed  procedure  is  as  follows: 

Step-1.  Pressure  gradient  across  each  branch  channel  for  the  given  geometry  without  orifice  is  obtained  by 
numerical  simulations. 

Step-2.  The  required  pressure  gradient  across  each  branch  channels  that  drives  the  desired  flow  is  calculated  using 
Hagen-Poiseuille  equation[18]  . 

Step-3.  For  the  pressure  difference  across  branch  channels  obtained  by  step-l&  Step-2,  the  orifice  diameters  are 
found  by  solving  the  orifice  equation  using  the  MATLAB  code. 

Orifice  equation: 

q  C.jaJIKp 

fp{  1~P) 

Where  p=d/D 

Step-4.  Numerical  simulations  are  performed  and  flow  deviation  is  measured  with  reference  to  the  desired  flow 
for  each  channel.  Based  on  the  percentage  deviation  from  the  desired  flow  across  each  branch  channel,  orifice  diameters 
are  determined  for  the  next  iteration. 

If  the  obtained  flow  rate  is  less  than  the  desired  flow,  orifice  diameter  is  increased  by  the  calculated  percentage 
deviation  and  vice  versa.  The  upper  limit  for  increasing  the  orifice  diameter  is  the  diameter  of  channel  itself  i.e.  maximum 
flow  possible  in  a  channel  is  without  any  orifice.  Detailed  procedure  of  arriving  at  the  pattern  flow  distribution  is  described 
in  the  flow  chart  as  shown  in  Figure  4. 

START  j 

Assign  required  volume  f low  rate  in  ever^ 
individual  channel  (Q  channej)  and  ploT  the 
required  f  low  pattern 

Based  on  the  Q  channei,  calculate  required 
pressure  drop  (AP  required)  in  every 
individual  channel  using  Hagen-Poi seuille  Eqn 


For  every  channel  Calculate  AP  =AP  actual  -  AP  reguired 


For  all  positive  values  of  AP  f  ind  f  irst  set  of  orif  ice  diameters 
using  MATLAB  code  for  respective  channels  at  reguired  f low 

Solve  Trial  1  case  (with  orif  ices)  and  obtain  pressure  drop 
across  everv  channel 

Compare  the  channel  pressure  drop  obtained  in  Trial  1 
case,  wiTh  AP  required  and  Then  f  ind  percenTage  deviation 
w.  r.  to  channel  pressure  drop  obtained  in  Trial  1 

FurTher  improvise  The  orifice  diameTers  obTained  in  Trial 

1  case  by  the  percenTage  deviation  obtained  in  previous 
steD  and  obtain  second  set  of  orif  ice  diameters 

Implement  the  second  set  of  orif  ice  diameter 
obtained  in  previous  step  and  solve  Xrial  2  case 


When  the  %  deviation  with  respect  To  expected  f  low 
reaches  optimum  value  beyond  which  no  significant 
improvernent  is  observed,  the  iteration  stops 

Q  STOP  j 

Figure  4:  Methodology. 


Solve  Natural  flow  distribution 
case  (without  orifices).  Based 
on  this  case  determine  actual 
pressure  drop  (AP  actual)  in 
every  individual  channel 

■ 
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4.  NUMERICAL  INVESTIGATIONS 

Figure  5  shows  the  computational  model  used  for  CFD  analysis.  It  consists  of  ten  number  of  channels  with  a  pitch  of 
70mm  for  the  flow  pattern-1  and  pattern-2.  First  9  channels  are  considered  for  the  flow  pattern-3.  Z-type  manifold 
configuration  is  considered  for  the  present  study.  Diameters  of  header  and  channel  are  taken  as  12mm  and  6.5mm 
respectively.  Total  length  of  the  manifold  considered  is  700mm. 

Water  is  considered  as  coolant.  Density  of  water  is  taken  as  998kg/m3for  the  numerical  simulations.  Bottom  inlet 
and  top  outlet  are  considered  for  the  manifold  system  and  3Dnumerical  simulations  are  carried  out  in  ANSYS  Fluent 
V17.0  CFD  solver.  Laminar  flow  with  yelocity  inlet  and  pressure  outlet  boundary  conditions  are  imposed.  Using  the 
automatic  grid-generating  tool  ANSYS  ICEM-CFD,  hybrid  grid  is  generated  for  the  whole  computational  domain.  Very 
fine  grids  are  used  in  the  some  areas  where  the  variation  in  the  flow  parameter  is  high,  for  accurate  prediction  of  fluid 
properties,  for  example  at  the  combining  and  dividing  junctions  of  the  manifolds  with  the  branch  channels,  element  sizing 
is  made  very  fine  (0.002mm)  to  accurately  capture  the  minor  losses  during  turning  of  the  fluid.  Average  element  quality  of 
the  entire  fluid  domain  is  maintained  at  0.7  for  all  the  analysed  cases.  The  analysis  was  performed  using  Steady  state  solver 
with  incompressible  flow.  The  fluid  flow  conditions  in  the  domain  are  governed  by  the  conservation  laws  of  mass  and 
momentum.  These  governing  equations  are  solved  using  pressure-based  solver.  Second  and  higher  order  numerical 
schemes  were  used  for  discretization  of  governing  equations.  The  SIMPLE  algorithm  is  selected  for  pressure-velocity 
coupling.  The  convergence  criterion  is  set  in  such  a  way  that  the  residual  of  each  calculated  parameter  is  less  than  10-5.  All 
the  simulations  were  carried  out  using  parallel  processing  system  having  32  core  Intel®  Xeon®  processor.  All  the  cases 
reached  the  converged  solution  within  about  1500  iterations  on  a  64  GB  RAM  work  station. 


Grid  independence  study  is  carried  out  for  the  cases  of  natural  flow  configuration,  that  is  the  configuration  with  no 
orifices  in  the  branch  channels.  Flow  patterns  are  compared  for  0.5,  0.85  and  1.24  million  elements  and  variation  in  the 
flow  and  pressure  drop  parameters  is  found  to  be  within  1-5%  as  shown  in  Figure  6.  Since  this  variation  is  small,  all  the 
cases  are  solved  for  the  grid  size  between  0.8  to  1.3  million  elements. 

5.  NUMERICAL  RESULTS  AND  DISCUSSION 

Flow  rates  across  each  branch  channels  for  all  three  patterns  obtained  by  numerical  simulations  are  presented  in  the  Figure 
7a,  Figure  7b  and  Figure7c.  Average  percentage  deviation  vs.  the  iteration  number  is  plotted  in  Figure  8  and  it  is  observed 
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that  no  signihcant  change  is  seen  from  second  to  third  iteration  for  all  the  three  patterns.  Hence  the  iterations  are  not 
continued  after  the  third  one.  Orifice  diameters  obtained  from  the  simulations  are  presented  in  Table-1.  It  is  to  be  noted 
that,  where  the  orifice  sizes  are  not  specified,  the  pipe  will  be  the  straight  channel  without  orifices. 


*  r.  □  t l  r □  I  f Idw  *  Requirtd  fKaw 


Channel  No. 


Figure  7a:  Pattern-1.  Figure  7b:  Pattern-2. 

- o Nn  turHltlDW  - A R±quTred  JTaw 


Figure  7c:  Pattern-3. 


Iteration  No 


Figure  8:  Average%  deyiation  from  there  quired  flow  vs.  Iteration  No. 


Table  1:  Orifice  Diameters  Obtained  after  Second  Iteration 


Channel  No. 

Pattern-1 

Pattern-2 

Pattern-3 

1 

2.4 

4.5 

1.6 

2 

- 

4.5 

1.6 

3 

2.3 

- 

- 

4 

- 

2.7 

1.6 

5 

2.3 

4.5 

1.6 

6 

- 

4.5 

1.6 

7 

2.2 

- 

- 

8 

- 

2.4 

1.5 

9 

2.1 

3.8 

1.5 

10 

5.7 

3.3 

- 
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6.  EXPERIMENTAL  SET-UP 

To  validate  the  analysis  results,  experimentation  is  carried  out.  Figure  9  shows  experimental  set-up.  For  a  better 
understanding,  detailed  schematic  layout  that  represents  all  the  components  of  the  experimental  set  up  is  shown  in 
FigurelO.  Inner  diameter  of  the  rota  meter  is  6.5mm  which  is  also  the  inner  diameter  of  the  channel.  Total  length  of  the 
manifold  is  700mm  and  height  of  the  channel  is  440mm.  Manifold  tubes  are  of  inner  diameter  12mm.  All  the  joints  are 
connected  with  the  leak  proof  quick  dis-connect  couplings.  Each  channel  rotameter  is  graduated  between  0-0.3  lpm  with  a 
least  count  of  O.Ollpm.  Accuracy  of  the  rotameters  is  ±0.  25%.  At  inlet  side  two  rota-meters  are  connected.  The  first  one 
has  range  of  0-5  lpm  with  least  count  of  0.2  lpm.  The  second  one  has  range  of  0-251pm  with  least  count  of  0.5  lpm.  A 
submersible  water  pump  of  201pm  pumping  capacity  connected  to  0.5kW  motor  that  runs  on  single  phase  power  supply  is 
connected  to  the  water  reservoir.  Water  reservoir  size  is  0.25mx0.25mx0.5m  with  30  litres  capacity.  Water  level  indicator 
is  attached  to  this  reservoir.  The  water  delivery  line  that  runs  from  the  pump  has  a  flow  control  valve  that  can  be  operated 
manually.  Since  the  motor  is  of  constant  RPM,  when  the  flow  control  valve  is  operated,  water  which  is  excess  to  the  test 
requirements  is  made  to  return  to  the  reservoir  through  a  bypass  pipe  line. 


Figure  9:  Experimental  Set-up  Figure  10:  Schematic  Layout  of  the  Experimental  Set-up 


Before  experimentation,  each  rotameter  is  verified  with  known  volume  flow  rate  of  water  to  find  the  correction 
factor.  For  each  pattern  flow  configuration,  experiments  are  conducted  for  4  times  to  verify  repeatability. 

7.  RESULTS  AND  DISCUSSIONS 

Figure  1 1  provides  a  comparison  of  flow  rate  through  branch  channels  between  the  results  obtained  from  simulations  and 
experiment.  It  is  observed  that  there  is  an  average  deviation  of  9.8%  between  the  two. 


Figure  11:  Results  Comparison  for  Natural 
flow-Experimental  &  Numerical 
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Comparison  between  the  required  flow,  results  obtained  from  CFD  simulations  and  experimental  results  are 
plotted  in  Figure  12a,  Figure  12b  and  Figure  12c  for  pattern-1,  pattern-2  and  pattern-3  respectiyely. 


Average  deviation  from  the  flow  requirements  and  the  experimental  results  for  pattern-1,  pattern-2  and  pattern-3 
are  found  to  be  4.5%,  3.9%  and  6.7%  respectively  which  are  in  close  agreement  to  each  other. 

8.  CONCLUSIONS 

In  the  emerging  new  requirements  of  cooling  of  the  electronic  units,  the  uniform  flow  distribution  may  not  always  provide 
optimal  solution.  For  the  similar  geometry  configurations  and  equal  spacing  of  the  electronics  inside  the  racks,  three  types 
of  patterns  are  taken  as  case  studies  and  solved  with  the  proposed  methodology.  The  methodology  developed  to  arrive  at 
desired  pattern  flows  is  demonstrated  using  numerical  simulations  and  the  same  is  validated  by  experimental  results.  The 
average  percentage  deviations  between  the  required  flow,  CFD  and  experimental  results  are  summarised  and  presented  in 
Table-2.  It  is  seen  that  all  are  in  close  agreement  with  each  other. 


Table  2:  Ayerage  Percentage  Deyiations 


Pattern 

Required  Flow 
&  CFD  Result 

Required  Flow  & 
Experimental 
Result 

CFD  &  Experimental 
Result 

1 

4.5 

9.16 

11.4 

2 

3.9 

9.3 

6.7 

3 

6.7 

4.4 

9.3 

This  methodology  provides  an  effective  and  simplified  solution  that  can  be  adopted  in  electronics  cooling  in  case 
of  dissimilar  heat  load  characteristics  with  repeated  patterns. 
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